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Abstract 

The Schrodinger and Heisenberg evolution operators are represented in phase 
space r*R" by their Weyl symbols. Their semiclassical approximations are con- 
structed in the short and long time regimes. For both evolution problems, the 
WKB representation is purely geometrical: the amplitudes are functions of a Pois- 
son bracket and the phase is the symplectic area of a region in T*R" bounded 
by trajectories and chords. A unified approach to the Schrodinger and Heisen- 
berg semiclassical evolutions is developed by introducing an extended phase space 
X2 = T*(r*M"). In this setting Maslov's pseudo differential operator version of 
WKB analysis applies and represents these two problems via a common higher 
dimensional Schrodinger evolution, but with different extended Hamiltonians. The 
evolution of a Lagrangian manifold in X2 , defined by initial data, controls the phase, 
amplitude and caustic behavior. The symplectic area phases arise as a solution of 
a boundary condition problem in X2- Various applications and examples are con- 
sidered. Physically important observables generally have symbols that are free of 
rapidly oscillating phases. The semiclassical Heisenberg evolution in this context 
has been traditionally realized as an h power series expansion whose leading term 
is classical transport. The extended Heisenberg Hamiltonian has a reflection sym- 
metry that ensures this behavior. If the WKB initial phase is zero, it remains so 
for all time, and semiclassical dynamics reduces to classical flow plus finite order h 
corrections. 



1 Introduction 



The version of quantum mechanics most suitable for semiclassical dynamics is Moyal's 
|l|] representation. This classical like statement of quantum mechanics is achieved by 
the Wigner-Weyl 0, |^ correspondence, which uniquely associates Hilbert space opera- 
tors with Weyl symbols (functions) on phase space. This correspondence transforms the 
non-commutative product of operators into a non-commutative * product of symbols. 
Expectation values are realized as integrals over phase space. 

For systems whose coordinate manifold is flat, the Wigner-Weyl isomorphism is most 
compactly described 0, |, Q in terms of the quantizer. Let T(x) = exp[z(j9-g — g-p)/^] be 
the Heisenberg translation operator. Here x = {q,p) are the position, momentum coordi- 
nates of the linear phase space xi = T*W^, and x = {q,p) are the corresponding quantum 
counterparts. The parity operator on Hilbert space is Pipiq) = ,4' ^ -^^(^g)- The 

product of these two unitary operators defines the quantizer, A{x) = 2"T(2x)P. In terms 
of A (a;), the one-to-one Weyl symbol-operator pairing is 

A = h-'' f dxA{x)A{x) , A{x) = [A]^ix) = TtAA{x) , 

where h = {2'Kh) . The first statement above is Weyl quantization; while the second is 
de-quantization. Whenever the operator is a density matrix p, the quantity /i~"p(a;) is 
the Wigner distribution. The above map, from A to A{x), gives the same result as the 
Wigner transform ( |A.3| ) of the kernel (gl^lg'). 

The star product is defined by the requirement that [AB = A* B cf. Appendix A. 
For small h, the * multiplication is approximately the commutative product of symbols, 

A * B{x) = A{x)B{x) + ih/2{A, B}{x) + 0{h^) . 

The Poisson bracket term measures, to leading h order, the non- commutative character 
of the star product. 

The two basic statements of quantum dynamics are the Schrodinger and Heisenberg 
pictures. On the space, L^(Mg), each self-adjoint Hamiltonian H generates a unitary 
Schrodinger evolution, U{t) = exp[—itH /h]. An initial density matrix po has Heisenberg 
evolution, p{t) = U{t)poU{ty. Let H, U{t) and p(t) be the Weyl symbols of U(t) and 
pit), respectively. In this setting the equations of motion are 

ihdU{t,x)/dt = H{x)*U{t,x), (LI) 
dp{t,x)/dt ={H,p{t)}M{x), (1.2) 



with initial conditions U{0,x) = 1, and p(0,x) = po{x). In ( |L2|) , the bracket operation is 
defined by 

{A, B}m = \[A,B]„ = \{A*B-B*A). 
in in 

This is the Moyal bracket. Like the quantum commutator, it is bilinear, skew and obeys 
the Jacobi identity. 
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It is evident that the Weyl-Heisenberg symbol evolution problem ( |1.2| ) is formally as 
close as it can be to classical dynamics. Since {A,B}m = {A,B} + 0{h^) , Bohr's cor- 
respondence principle is realized in a transparent form. The time dependent expectation 
value of an observable A is conveniently given by the phase space integral 

{A)t = Tr Ap{t) = h-"" I dx A{x)p{t, x) . (1.3) 

The main goal of this paper is to obtain both short and long time WKB approximate 
solutions for the Weyl symbol evolution problems ( |1.1| ) and (|1.2|) . In Theorems 1 and 



2, we prove that the phases entering these semi-classical representations are invariant 
geometrical quantities — namely the symplectic areas defined by certain closed loops in 
phase space. Furthermore, the amplitude functions are given by determinants of a Poisson 
bracket. 

We employ a common method to solve both problems (|1.1|) and (|1.2| ). Standard 



Weyl symbol calculus identities show that H* and {H, ■}m are pseudodifferential opera- 
tors, in other words they are functions of V^: and multiplication by x. If one introduces 
an extended phase space X2 = T*{T*W^), then Maslov's ^ WKB analysis applies. In 
particular, the semiclassical approximations for the evolutions U (t, x) and p{t, x) are char- 
acterized by a Lagrangian manifold, and the time evolution of this manifold is determined 
by an extended Hamiltonian in X2- In order to distinguish these two phase spaces we call 
Xi the primary phase space (PPS), and X2 the secondary phase space (SPS). Projecting 
the X2 analysis back onto the primary phase space provides explicit WKB representations 
of U{t,x) and p(t,x). 

The results in the literature nearest to this paper are found in three seminal works: by 
Karasev and Nazaikinskii p|. Berry |^, and Marinov JlOl. In the first of these, an analyt- 



ical approach to a generalized semiclassical expansion based on a SPS with its associated 
left-right representation of *-product and Hamilton- Jacobi equation was developed. In 
the second pair of works, the geometrical (symplectic phase) interpretation of the Wigner 
eigenfunction and of the dynamical WKB approximation was achieved. Specifically, in 
the small time sector, where the time evolving X2 Lagrangian manifold remains single 
sheeted, Marinov found the WKB amphtude and phase for U(t,x). In addition to pro- 
viding an entirely different proof, our treatment extends these known results in several 
ways. First, we determine the long time version of the WKB approximation where the X2 
Lagrangian manifold may be multi-sheeted with several points corresponding to a given 
X e Xi- Second, by incorporating Maslov's results one maintains control of the error 
estimates. Third, we solve a generalized version of ( |1 . 1|) with rapidly oscillating initial 
data. This latter generalization arises for the evolution U(t + to), to > 0. At t = the 
initial data generically will have the rapidly oscillating form, A'"o(x)e**''*^^^/'*. 

In spite of the importance of problem ( |1.2D , it has received insufficient study. The 



version of the Heisenberg problem that has been investigated in detail [11, |1^, 0, n 
occurs for the case where the observable A is semiclassically admissible, namely its Weyl 
symbol is h dependent and allows a finite order small h asymptotic expansion cf. (p.6|). 



In |T^ a connected graph method was constructed to obtain the coefficients of the h 
asymptotic expansion of [U{tyAU{t)]vj{x). Based on this latter formalism, numerical 
calculations have shown that the /i^-order semiclassical expansion accurately describes 
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noble gas atom-atom scattering. However, semiclassically admissible symbols are not 
suitable for representing [|1^ density matrices. The requirement that one has a pure 
state, Po * Po = Po, means that po{x) is an /I ^ rapidly oscillating function cf. ( p.l2| ). 
The paper has the following organization. Section 2 presents the extended Schrodinger 



equation in the SPS setting and shows how the problems (|1.1|) and (|1.2| ) are particular 
examples of this extended evolution. Also, the geometric, commutative and coordinate 
relationships between the PPS and SPS are reviewed. The trajectories entering the WKB 
expansion are determined by the solutions of a boundary condition problem. In Section 3, 
a detailed analysis of the BC problem and the manner in which it constructs phase space 
loops is carried out. The semiclassical expansions for U{t,x) and p{t,x) are consolidated 
in Theorems 1 and 2. Finally, in Section 4 the group generated phase addition rules are 
obtained; the mutual compatibility of the Schrodinger and Heisenberg semiclassical evo- 
lution is established; the exact solutions for quadratic Hamiltonian systems are compared; 
and, the ^ — > asymptotics of [U{t)AU(ty]^{x) for semiclassically admissible operators 
A is realized as a special case of Theorem 2. There are three appendices containing spe- 
cialized aspects of Weyl symbol calculus, pseudo differential operators, Jacobi fields, and 
a Poincare-Cartan identity. 



2 Secondary Phase Space Dynamics 

In the first part of this section we show how to interpret the H* structure in the equations 
of motion as the action of a pseudodifferential operator (\EfDO). The normal ordered 
symbols of these operators are then functions on the secondary phase space. For the 
Schrodinger and Heisenberg problem we show how the classical flow in X2 is composed in 
terms of the standard Hamiltonian mechanics on xi- The necessity of representing the 
initial value form of the density matrix as a rapidly oscillating h function is confirmed by 
constructing the semiclassical Weyl symbol for a pure state. 

Introduce the SPS canonical operators X, Y acting on functions f{x): Xf{x) = xf{x) 

2 1 

and Yf{x) = —iKS7xf{x). Let 7i(X, Y) be a normal ordered, symmetric \E'DO defined by 
the symbol T-C{x,y). On the coordinate domain M^"^, the extended Schrodinger evolution 
problem is 

ihd^{t,x)/dt = H{X,Y)'^{t,x) . (2.1) 

The two problems ( |1 . 1| ) and ( |1.2| ) are particular cases of (|2.1| ). The appropriate realization 
of TC in each case is determined by employing the left-right * product representation ||, 
(see Appendix A). Setting H{x,y) = Hi = H{x — \ Jy) selects the Schrodinger problem 
( pTTl) , whereas the choice 0,2 = H{x — \ Jy) — H{x + \ Jy) gives the Heisenberg problem 
(|1.2|). The solution \E'(t, x) is either U(t,x) or p(t,x) depending on the choice of H. In 
SPS, the scalar function ?i is a Hamiltonian which is the generator of X2 classical flow. 

Consistency of this extended Schrodinger problem requires that the normal ordered 
symbol H{x ± ^Jy) define a symmetric operator on L^(M^"'). That is does so is a con- 
sequence of the argument structure (x ± \Jy) which implies that cf. ( [A.16[ ) the normal 



ordered and Weyl ordered symbols are the same function. Since the Weyl ordered symbol 
is real, one has that the ^DO Hamiltonians in ( |2.1D are formally self-adjoint. 
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Consider the coordinate systems of xi cind X2 and their interconnections. Let y be 
the momentum variable in SPS so that a general point m G X2 has Cartesian coordinates 
z{m) = {x,y). The symplectic structure of Xi and X2 is determined by their respective 
canonical 2-forms 



id 



= dp A dq = \ Jai3 dxp A dxa , J 



/„ 

-In 



= dy A dx = \ Jjk dzk A dzj , 

where In is the n x n identity matrix, and J is the symplectic matrix having the block 
form of J with substituted for /„. Given the 2-forms, cu^^^ and a;^^-*, the corresponding 
Poisson brackets are 

{f,gUx) = Vf{x) ■ JVg{x) , {F.GUz) = VF{z) ■ JVG{z) . 

The coordinate functions, x = {q,p) and z = {x,y), are canonical in xi and X2, namely 
{xa,x^}i = Jap and {zj,Zk}2 = Jjk- In order that be a standard Schrodinger 

equation the coordinates x = {q,p) must commute. This is the case since {xa,Xf3}2 = 0. 

As one sees from functional forms of 7^, a second natural coordinate system for X2 is 
defined by the (/,r) variables 

l = x-\jy, r = x + \Jy, /, r G M^" , 

(2.2) 

x = ii^{l + r), y = J{l-r). 

The link between the coordinate systems (x, y) and (/,t), and their commutative struc- 
tures is the following. 

Lemma 1. (i) The SPS commutative properties of the projections l,r : M^" —>■ M^" are 

{J'aJl3}2 = Ja/3 , a, p] 2 = — J aP ■, {la-,'^ p}2 = ^ ■ 

(ii) The I variables are Poisson; the r variables are anti-Poisson, namely 

{f ol,go l}2 = {/, g}i o /, {f or,go rja = -{/, g}i o r. 

(iii) The X2 Poisson bracket and symplectic form have the I x r decomposition 

{F,G}2 = {ViF)^J^p{ViG)p - {VrF)^J^p{VrG)p, 
Jjk dzk A dzj = Jap dip A dla — Jap drp A dva ■ 

The space structure for the left-right variables is xi ® Xi 3 Here Xi denotes the 

phase space defined by the 2-form —\Jap drpAdva. Often the space Xi is labelled 
the double phase space. We use the name secondary phase space for cotangent bundle 
X2 = T*{T*R^) in order to distinguish it from xi ® Xi- The left-right map, ( |2.2|) is 
diffeomorphic. 

The basic ingredients in any semiclassical approximation are the classical action phase 
and the amplitude functions determined by classical transport. Because the extended 
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Hamiltonians are the simple composite functions H o / and H o I — H o r, it is possible 
to describe the X2 phase space motion as appropriately weighted sums of Xi trajectories. 
Let g{T\xo) = {q{T\xQ); p{t\xq)) denote the solution of the PPS Hamiltonian system 

g{r\xo) = JVH{g{r\x^)) , g{Q\xo) = Xo . (2.3) 

The SPS motion is labelled G{t\zq) = {x{t\zo); y{T\zo)) , where 

G{t\zo) = JVn{G{T\zo)) , G{0\zo) = z, . (2.4) 



The initial data zq for system ( |2.4| ) also has the left-right representation Iq = l{zo) ,ro = 
r{zo). 

Lemma 2. The solution of the SPS Hamiltonian system, in terms of xi flows, takes the 
form: (i) For the Schrddinger problem with Tii = H o I 

x{t\zo) = ^{giT\lo) +ro) , yir\zo) = J{g{T\lQ) - ro) . (2.5) 

(ii) For the Heisenberg problem with 7-^2 = Hol — Hot 

^i'^lzo) = ^{girllo) + giT\ro)) , y{T\zo) = j{g{T\lo) - g{T\ro)) . (2.6) 

Proof. Consider part (i). Using the Poisson equations of motion, it is convenient to state 
the SPS flow in the I, r variables. Lemma 0(iii) implies 

Ut) = Hol}^ = J^^VfsH{l{T)) , r,(r) = {r„, Hol}, = 0. (2.7) 



So 1{t\Iq) = 5'(r|Zo) and the r-motion is a constant, r(r|ro) = tq. Using ( |2.2| ) recovers 
( p.5| ). A similar argument verifies part (ii). □ 

The identities ( p.7| ) show that the left-right motions are decoupled. This decoupling 
occurs because {L(t), T/3(r)}2 = 0. Further note that the system Hi = H o I has 2n 
constants of motion; ra{T) = const for a = 1 . . . 2n. Nevertheless, Tii is not completely 
integrable since the functions r^^T) are not in involution, i.e. {?"a(T), r^(r)}2 = —Ja/3 7^ 0. 

Lagrangian manifolds in SPS create canonical transformations in PPS. Let us clarify 
how this is realized in terms the left right projections. Suppose A is a X2 manifold which 
is locally defined by the smooth phase S : U C R^" — * M., i.e. y = VS{x),x G U. This 
identity relates / to r by 

r-/ = JV^(l(/ + r)). (2.8) 

If det(/ ± \ JS") , X E U then the implicit function theorem defines I = Z(r) and 
r = r{l). These transformations are canonical. This follows from 

dl (/ - Us") 



dr {I + \JS")' 

The right side above is the Cayley transform of the symmetric S" and so defines a sym- 
plectic matrix. A systematic study of generating functions of this type is found in 
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An important aspect of WKB analysis is the h singularity structure of the Cauchy 
initial data. As noted in the Introduction, a pure state density matrix must be a rapidly 
oscillating function. Let us clarify this situation with an example. Suppose ipolq) = 
n{q) exp[is{q)/h] is a unit normalized L^(M^) wave function, which defines a rank one 
density matrix, po = |V'o)(^o|- The functions n, s are real. The Weyl symbol of po is given 
by the Wigner transform ( [A.3|) 



Po{x;h) = / dvn{q + ^v)n{q — ^v) exp ^(j){v,x) , (2.9) 

(p{v,x) = s(g + |f ) — s(g — |f ) — p ■ f . (2-10) 

The ^ — > asymptotic form of po may be calculated by a stationary phase approxi- 
mation [0. Let (j)'{v,x) = Vv(f){v,x) and (f)"{v,x) = V^V^0(f,a;). The critical points of 
satisfy 

2p = Vs{q-lv) + Vs{q+lv). (2.11) 

Suppose {vo,Xo) is a solution set for this equation where det ^"(fo, a;o) 7^ 0, then the 
implicit function theorem defines a function v = v{q,p) = v{x) obeying ( p.ll| ). Since 
4>'{v, x) is even in v, the roots occur in pairs. If f = v{x) is a solution, then so is w = —v{x). 
Incorporate this pairing behavior via the notation: q± = q ± v{x)/2, p± = Vs(g±), 
x± = [q±,p±). In the non-singular region where det (f)"{v{x),x) 7^ 0, the Wigner transform 
above has asymptotic form 

po{x]h) ^h^\det 0"(i;(x),x)|-^/2n(g+)n(g_)2cos|^5o(x) + ^sgn</)"(t;(x),x)| . (2.12) 

The phase 5*0(0;) = (j){v{x),x) and sgn denotes the signature of a symmetric matrix. If x 
is in the neighborhood where ( |2.11| ) has no solutions, then po{x) is 0{h°°) small. 

As Berry first noticed [^] (in the n = 1 version of this problem), the phase 5*0 is a 
symplectic area. To see this in the present context, define the Xi Lagrangian manifold 
A = {x G Xi\p = Vs(g)}. Construct two paths between the end points x_ and x+. Let the 
first path, C{x-,x+), be a straight line from x_ to x+. As the second path, 7a(x+,x_), 
choose any curve from x+ to X- lying on the surface A. Then one finds 



So{x) = (t p-dq. (2.13) 



The geometry of the loop in ( p.l3| ) is the well known chord mid-point construct |]T0 



2l| , p2| , The end points x± of the directed chord C(x_,x+) lie on A, and x is the 
mid-point of this chord. By Stokes theorem, So{x) is also the area of any membrane 
having boundary C{x-,x+) + ■yx{x+,X-). 



3 Symplectic Areas and WKB Phases 

This section shows how the secondary phase space WKB approximations for U{t, x) or 
p{t, x) can be reformulated in terms of primary phase space flows. This reduction process 
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transforms all the Xi WKB phases into symplectic areas defined by closed loops in x\- 
A key element of this analysis is the solvability of the appropriate two point boundary 
condition (BC) problem. We pose the EC problem in terms of Lagrangian manifolds and 
construct both short and long time solutions. 

To begin, we summarize Maslov's WKB expansion for the extended Schrodinger 
problem (|2.1| ). The t = state for system is assumed to have the generic form, 

\E'o(x) = ao(3;)e***-^^/^. The real functions ao,/?o have support on the domain C M^"-. 
The phase /3o defines a X2 Lagrangian manifold Aq, via y{x) = V/5o(x),x G -Do- Let 
Ili{x,y) = a: be the projection onto M^". It is assumed that projection Hi : Aq — Dq is a 
diffeomorphism . 

To proceed three assumptions are required: (a) The function 7i G T^(M^"),m > 2. 
cf. Appendix A. (b) For arbitrary finite time intervals [— T, T] , the SPS Hamiltonian 
system (|2.4|) has unique smooth solutions, i.e. G(t\x,y) G C°°{[—T,T],'Rf^). (c) Let x 
be a non-focal point {cf. Definition 3) with respect to the Hi projection of the manifold 
A( = G(t)Ao. Assume there are a finite number of points ocq G -Dq, j = 1, . . . , N such that 
U,Git\xlVl3oixi))=x. 

2 1 

Each of these assumptions is required for an evident reason: (a) ensures that Ti.{X, Y) 
is a well defined \l/DO, (b) prohibits finite time runaway trajectories, and (c) assumes the 
existence of one or more solutions to the BC problem: which trajectories starting from 
Aq have x as their final coordinate position? Locally varying {t, x) in (c) determines the 
initial xq as a function of {t, x), namely Xq = xl{t, x). 

Under the hypotheses (a-c) the h —>■ asymptotic solution of the Cauchy problem 



(p. ID, with initial data \E'o, is 



^{t,x) = ^'''{t,x)[l + 0{h)], (3.1) 
vE'-(t,x)= f;-Ml^exp{^/5,(t,x)-ym,(t,x)}, te[-T,T]. (3.2) 

Here mj{t,x) is the Maslov index for trajectory {G{t\zq) : r G [0,t]}, where Zq = 
{xl{t,x),'VPo{xl(t,x))) G Aq. The phase and amplitude functions are 

f3,{t, x) = (3o{xi) + r [y{T\zi) ■ x{T\zi) - H{G{T\zi))] dr , (3.3) 







V dxl J 

0,(t,x) = (^exp j^trV,Vyn{G{T\zi))dT^ao{xi). (3.5) 



2 1 



The family of operators T-C{X,Y) consistent with the class T™(M^") is large. The 
growth index, m, must be 2 or greater so as to include kinetic energy in H. The operators 

2 1 

?i(X, y) need not have a polynomial dependence in momentum Y, nor do they need to be 
partial differential operators. This generality is essential to the SPS method we employ: 
normally, the \E'DO defined by Tii and 7-^2 will not be a finite order partial differential 
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operator. The phases l3j{t,x) have a geometric meaning. Define the Ti. fiow transported 
Lagrangian manifold, Aj = G(t) Aq. Then the surfaces yj{t,x) = V(3j{t,x) describe the 
different sheets of A^. Representation (|3.lH O|) is Theorem 10.5 of 0. 



In phase space quantum mechanics, the range of physical systems one may describe is 
controlled by the allowed functional form of the Weyl symbol Hamiltonian. For reasons 
of notational simplicity we presume H is static. Including time dependent Hamiltonians, 
as a subsequent generalization, is a straightforward matter. The WKB summary above 
indicates that two restrictions on H are necessary. 

Assumption Al. The Hamiltonian operator H is semiclassically admissible, specifi- 
cally H G T™(M^"') ,m > 2. A2. On finite time intervals [— T, T], H generates unique 
classical fiow, g{t\x) G C°°([-T, T], Mf). 

Assumption A is common to all the Propositions and Theorems of this section and 
will not be explicitly cited in them. 

Since I and r are linear functions of Assumption Al implies that the composite 
functions Tii = H o I and 7^2 = H ol — H o r are in T™(M^"), m > 2. Likewise, Lemma 2 
shows that Assumption A2 guarantees that y) G C°°([— T, T], M^") for both Tii and 

7^2 fiow. The requirement that H G T^(M^'i) (c/. Appendix A) means that H = H{x; h) 
has the h ^ uniform asymptotic expansion 

J 

H{x; h) = Ho{x) + J2 ^^Hj{x) + 0{h-^^^) . (3.6) 

3=1 

For many physical systems H{x;h) = Hq{x). For example, this occurs for a mass m 
charged particle moving in an external electromagnetic field. The higher order terms 
Hj[x) have no effect on the leading order WKB approximation. So from here on we will 
identify H with Ho and omit the subscript. 



3.1 The Boundary Condition Problem 

The assumption A says nothing about the solvability of the EC problem. Rather than 
dealing with this by assumption, as was done in the representation (|3.1| - p.5|) , we will find 
existence proofs and explicit formulas for the EC solutions. 

We discuss the EC conditions for the evolutions U (t, x) and p(t, x) in tandem. For the 
U {t, x) evolution, take the initial state of the system to be 

Uo{x) = Aro(a;)e^*o(^)/'' , x G Z^o ^ M^" (3.7) 

where the domain, Dq, is a simply connected open set. The phase $o defines the single 
sheeted Lagrangian manifold Aq = {z E X2\y{x) = V^o{x) , x G -Do}- The initial state for 
density matrix evolution p{t, x) is presumed to be 

Po{x) = ao(x)e'^« , xeDqCR'''. (3.8) 

Again, Aq will denote the simply connected manifold generated by Sq. The phases $o and 
^0 are both assumed to be C°°(M2n)_ 
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BC Problem. Determine the trajectories G{t\xo,yo) which begin on Aq and have the 
final position x & Df = HiAf . 

(i) The Schrodinger case: Find the functions Xq = Xo(t,x) and yo = yo{t,x) from 

\ \9{t\xQ - iJyo) + (xo + I Jyo)] = X , (3.9) 
yo = V%{xo), xoeDo. (3.10) 

(ii) The Heisenberg case: Find the functions Xq = Xo{t,x) and yo = yo{t,x) from 

\ [^(^ko - i^2/o) + 5'(^ko + \ Jy^)\ = X , (3.11) 
l/o = VSo(xo), XoeDo. (3.12) 

Lemma 2 justifies writing the left sides of (|3.9|) and (|3.11|) as the weighted sum of 
the left and right flows. It is natural to interpret each of these equations as a midpoint 
condition. Conditions (|3.1CI|) and ( p.l2| ) ensure that (xo,yo) lies on Aq. Observe that 
the initial density matrix ( |3.8| ) is not real valued like approximation ( [2.121 ). However, by 
superposition with the complex conjugate of (|3.8|) , the initial state can be made real. 

We treat the small time and big time solutions to the BC problem with different 
methods and assumptions. The small time problem may be formulated as follows. Notice 
that the midpoint and Aq boundary conditions may be combined into a single statement, 
Mt{x') = X. 

Definition 1. (i) In the Schrodinger problem set 

Mt{x') = Mt{x') = \ [g{t\x' - i JV$o(x')) + {x + \JV%{x'))\ = x . (3.13) 

(ii) In the Heisenberg problem set 

Mtix') = Mt{x') = \ [g{t\x' - \JVSo{x')) + g{t\x' + \JVSo{x'))\ = x . (3.14) 

The initial phase f3o{x) in the extended Schrodinger problem is equal to either ^o{x) or 
5*0 (x) depending on whether one has case (i) or (ii). Here it is assumed that the support 
of (3o{x) is all of M^", and that x is any point in M^". If the non-linear function Ait is 
invertible, then the BC problem has a single trajectory specified by the initial values, 
xo{t,x) = M.t~^{x) and yo{t^x) = V/3o(xo(t, x)). The arguments of the flows g{t) are the 
left-right projections of Aq. 

Proposition 1. Let the Hessians H"{x) and (3q{x) have x-uniform bounds 

\\H"ix)\\<c,<oo, \\(3'^ix')\\<2, (3.15) 

then there exists ati > such that the midpoint map Ait '■ is a diffeomorphism 

fort E [— ti,ti]. Specifically, xo{t,x) = Ait^^{x) is a unique C^{[—ti,ti],M?^) solution of 
the BC problems (i) and (ii), respectively. 

Proof. Consider the Heisenberg problem. For parameters {t,x) G ([— T, T], R^"), define 
the family of maps T^,^ : M^" M^n^ 

TtJx') = x + x'-Mt(x'). 
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The idea of the proof is to show that Tt^x is a contraction mapping. Its unique fixed 
point, x), obeys x = Mt{x*{t,x)) and estabhshes that Mf is invertible. Connect an 
arbitrary pair of points {xi,X2) by the hnear path x{C,) = Xi + C,{x2 — Xi). Integrate 
dTt^x{x(i))/di on ^ G [0, 1] to show 



Jo 



Whenever the integral above is less than one for all {xi,X2), then Tf ^, is a contraction 
mapping. Taking the x' gradient of Tt^x{x') it follows that 



|VT,.(x')|| <l 



\I - VgmW + 11/ - Vg{t\r')\\ + ^||V^7(t|0 - Vg{t\r')\\\\S'^ix' 



< 2(exp ci| t| - 1) . 

To obtain the final inequality, ||>S'o(a;')|| was replaced by its bound 2, and the Lemma ^ 
(Appendix A) estimates of the Jacobi field Vg were employed. Let ti > satisfy 
2(exp (citi) — 1) < 1, then for times |t| < ti the rightmost combination above is less 
than one and Tt^x is a contraction mapping for {t,x) G [— x M^". An application of 
the implicit function theorem to Mt(xo(t,x)) = x establishes that Xq G C^{[—ti,ti],M?^). 
The same argument applies to the Schrodinger BC problem. □ 

The merit of the short time BC solution given in Proposition 1 is that it is global — 
applicable for all {t,x) G [— x M^", and so in this time interval A^ has a diffeomor- 
phic x-projection. Nevertheless, the results are restrictive. Typically only a small time 
interval [— ti,ti] is allowed, and the Hamiltonian and initial phase must have no more 
than quadratic growth as | x| — oo. 

The BC problem in the large time regime has a different geometry. Let Aq have an 
x-projection. Do, that is contained within some compact set of M^". As Aj evolves away 
from Aq, it may develop into a multi-sheeted manifold. Here, it is helpful to recall some 
of the terminology related to caustic behavior of Lagrangian manifolds, A C X2- 

Definition 2. A point m G A zs non-singular if it has a neighborhood that is diffeo- 
morphically projected onto M^". // the opposite holds, it is singular. The set S(A), of 
all singular points, is called the cycle of singularities. Its complement is the regular set, 
7^ = A/S(A). 



Definition 3. The x projection of the cycle of singularities is the caustic or focal set, 
ni(S(A)) = A^(A). Let Dt = IIiAj denote the x-image of At, and }C{At) the associated 
caustic set. The non-focal set is defined as the complement of the caustic set. Aft = 
A//C(Ai). 

If At has more than one sheet, the inverse of projection Hi : A^ Dt ceases to be single 
valued: there are x & N't with multiple pre-images, {m^}^ C A^, satisfying Ilim-' = x. 
The result below is stated for the Heisenberg BC problem, with obvious adjustment it 
applies to the Schrodinger problem. 
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Proposition 2. Let x E Mt ^ Dt he a non-focal point with pre-images {m^ G Af}^. 
Then the Heisenberg BC problem has solutions 

Mt{xiit,x)) = x, J = !,■■■, N. 

The functions, a;g(t, x), are smooth and have support in a neighborhood, Mt '^Uj 3 x. 

Proof. Let Ft{x',x) = Mt{x') — x, where Ft is smooth on the domain {x',x) ^ Dq x 
Aft- The implicit function theorem (IFT) solutions generated from Ft{x',x) = 0, will 
construct x^t, x). In order to apply the IFT one must have a solution pair {xq, x), where 
Ft{xo,x) = and det'Vx'Ft{xQ,x) ^ 0. Consider one pre-image of x. Each G 
has a 1-1 association with a point in Dq^ via Xi^ = niG(— tim-'). The set {xo,x) is a 
solution pair of Ft. It remains to show that the determinant is non-zero. Since x G A/j, 
Lemma |10| implies that detVM((a;o) 7^ 0. Furthermore Vx'Ft{x',x) = VMt{x') shows 
that detVx'Ft{xo,x) 7^ 0. The IFT then estabhshes the existence of XQ(t,x) defined on a 
neighborhood, Uj 3 x. If Ft G C^,p > 1, then Xq G C^. □ 

Once A( has become multi-sheeted, typically the regular set TZt is divided by the cycle 
of singularities into disjoint subsets {A-j.}^ . On A^ C TZt, the restricted projection IljA^ = 
Dl C N't is diffeomorphic. In this circumstance, oc^Q{t,x) has the explicit representation 

xi{t,x) = UiG{-t\{U{y\) , xeDi. (3.16) 



Identity ( p.l6|) determines x^it^x) directly without requiring an IFT construction. 



As is evident from J{x,t) = det VA^t(x), cf. (|3.4| ), the gradient structures in the BC 
problem determine the WKB amplitudes as functions of the left and right Jacobi fields, 
Vg{t\l') and Vg{t\r'). It is of interest to find a more geometric representation of these 
amplitudes. This is achieved by constructing a suitable X2 Poisson bracket related to the 
BC problem. First characterize the initial manifold Aq by the (independent) constraint 
functions 

C(x', y') = y'- V%{x') , 7^{x', y') = y' - VS^{x') . 

The conditions ({x',y') = and ri{x',y') = define Aq in our two problems. A second 
set of functions must necessarily carry information about the Ti dynamics. Select the 
midpoint flow functions x{t\x',y') for this role. Then one finds 

Lemma 3. In terms of the Aq constraint functions, ( and f], the Tii and 7^2 amplitudes 
have the Xi Poisson bracket representation 

detVMt(x') = det{x(t),C}2(a;', V<l>o(x')) . (3.17) 
detVMt(x') = det{x{t),r]}2{x',VSo{x')) . (3.18) 

Proof. Consider the 7^2 case, an elementary calculation shows that 

{xitUvMx',VSoix')) = l[\/gmiI-ljS'^ix')) + Vgit\r')iI+^JS'^ix'))] 



[VM,(x') 



\af3 ■ 
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The determinant of this matrix is ( p.l^ ). The same type of argument gives ( |3.17| ). □ 

This type of bracket representation of the amphtudes is common ||19|, ^ in the SPS 
context. The components of constraint C (and rj) are commutative with respect the X2 
bracket. This involution property is a necessary and sufficient condition [2^| for the 
manifold Aq to be Lagrangian. In addition, the midpoint variables x{t) are commutative. 
This is a result of the fact that the functions x(t) take their values on PPS plane which 
is X2 Lagrangian. 

Our boundary condition analysis is concluded by identifying a xi loop structure that 
is created by each BC solution. These loops are constructed from two types of Xi liiis 
segments — chords and trajectories. The quantity C{xi,X2) denotes the chord (the R^" 
geodesic) from xi to X2. Trajectory type line segments are specified by 



T^{xi,X2) 

T-{xi,X2) 



{9{t\xi] 
{9{r\xi) 



T : 




r 



t)}, X2 = g{t\xi), 
-t)}, X2 = gi-t\xi) 



The ± labelling distinguishes between forward and backward evolution. 

Consider the Schrodinger problem with mid-point BC solution xo{x,t). The Aq man- 
ifold requirement ( p.lO|) , when stated in the left-right variables, is 



lo{t, x) = xo{t, x) - iV$o(xo(t, x)) , ro(t, x) = xo{t, x) + iV$o(a;o(t, x)) . (3.19) 

After time displacement t, the Tii fiow moves Iq = lo{t,x) and tq = ro{t,x) into the new 
positions k = ltit,x) = git\lo) and rt = rt(t,x) = vq , cf. Lemma 2(i). 

The pair of points {lt,ro) define the chord C(/t,ro), having midpoint x. However, 
{It, ro) are also the endpoints of the linked line segments, C(ro, lo) + T~^{lo, It). Altogether, 
for each x E Aft, one has (see Fig. 1) a three sided oriented loop 



L{t, x) = C(ro, /o) + T+(/o, h) + C{lt, ro) . 




Figure 1. The Schrodinger loop. Figure 2. The Heisenberg BC solution. 

3.2 Evolution Representations 

Based on these solutions of the BC problem the WKB approximation for U {t, x) takes 
the following form. 
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Theorem 1. Suppose x is a non-focal point associated with the manifold At. Let {L^{t,x)}i 
be the family of loops generated by the Schrodinger BC solutions {a;o(t, x)}^ , x E Aft- Then 
the 0{h) approximation ofU{t,x) is 

N 

f/-(t,x) = 5^iV,(t,x)exp|^$,(t,x)-ym,(t,x)}, (3.20) 
^j{t,x) = $o(x^o(^>a;)) + (I p ■ dq - H{li{t,x))t , (3.21) 

JL3{t,x) 

Nj{t,x) = \ detV Mt{xi{t, x))\'^/'^ No{xi{t, x)) . (3.22) 

Before proceeding to the Heisenberg case, it is useful to make a number of remarks 
specific to Theorem 1. An important special case occurs in the problem where f/(0, x) = 1. 
Here $0 = and Aq C M^". In particular, Iq = xq = tq, so the loop U{t,x) simplifies 
into one chord C{lt,lo) with midpoint x and one dynamical curve T'*"(/o,/t) connecting 
the chord end points. If A^ has one sheet, then = 1, and the theorem above gives 



the representation obtained in [10|. However, it should be noted that Marinov did not 



investigate the BC problem. Instead his approach was to construct the phase $ = $i(t, x) 



by solving the Hamilton- Jacobi (H-J) equation associated |T0|, |, ^ with the evolution 
problem (|1.1|) , 

dt^{t, x) - H{x - I JV$(t, x)) = 0, (3.23) 

with the Cauchy initial data ^{0,x) = 0. Because of the highly non-linear character of 
H{x) this H-J problem is much more difficult to solve than the BC problem ( p.l3|) treated 
above. It is straight forward to verify that $(t, x) given by ( p.21| ) is a solution of ( p. 23 ). 

Marinov's work established the central role of the phase ^{t,x) in the semiclassical 
approximation of U{t,x). However, knowledge of the H-J equation in the form ( p.23|) . 



and the use of ^{t,x) as the generator of a canonical transformation, cf. ( |2.8| ), goes back 
to Poincare . 

Quantum operator identities induce corresponding Weyl symbol identities. A simple 
example of this is the unitarity property, U(t) = f/(— t)^. Its Weyl symbol image is 
U{t,x) = U{—t,x)*, which in turn requires that $(t, x) = — $(— t,x) and N{t,x) = 
N{-t,x). 

The generalization of Theorem 1 to non-static Hamiltonians H{t, x) is achieved if the 
phase ^j{t,x) is modified by the replacement 

H{li{t,x))t^ f H{g{T\li{t,x))dT. 
Jo 

Adjusting the notation above to cover the Heisenberg problem is easy. In this case 
lo(t,x) and ro(t,x) are defined by (|3.19| ) with 5*0 replacing $o- The right point rt is no 
longer the fixed point tq, but instead is rt = rt(t,x) = g(t\rQ(t,x)). The curve between 
and ro is T~(rt,ro). The resultant BC loop is the four segment structure (see Fig. 2) 

W{t, x) = C(ro, /o) + T+(/o, h) + C{lt, rt) + T-{rt, tq) . 



13 



The loop W(t,x) is like the Schrodinger BC loop L{t,x), but has an additional side 
T^{rt,ro) inserted at the point r^. However, there is another simpler loop generated by 
the BC solution. Denote by 7(t|C(ro, /q)) = {^'(^l^^') • ^' ^ C'(ro,/o)}- This curve is the 
flow translated image of the chord C(ro, Iq). In general it is no longer a straight line but 
it has end points (rtJt), so 

L{t,x) = C{k,n) + j{t\C{ro,lo)). 
defines a two segment closed curve, see Fig. 3. We refer to L(t, x) as the Heisenberg loop. 





Figure 3. Heisenberg loop area L(t,x). Figure 4. S{t,x) as an evolving front. 

Theorem 2. Suppose x is a non-focal point associated with the manifold At. Let {L^it, x)}^ 
be the family of Heisenberg loops generated by the BC solutions {xQ{t,x)}^ ,x G N't- Then 
the 0{h) approximation of p{t,x) is 



N 

p'''{t,x) = ^aj(t,x)exp |^S'j(t,x) - ymj(t,x)| 



Sj{t,x) =So{xl{t,x))+ <h p-dq, 

J L3{t,x) 

aj{t,x) = \ detV Mt{xi{t, x))\'^/^aoixi{t, x)) . 



(3.24) 

(3.25) 
(3.26) 



Proof. Within the present setup, direct computation establishes both Theorems 1 and 2. 
Consider the Heisenberg case. The dr integral contribution to the X2 WKB phase in (|3.3|) 
has both a y ■ x and a 7^2 contribution. The first of these is 



{9{t\Io) + 9{T\ro)) ■ J{g{T\lo)-g{T\ro))dT (3.27) 



p ■ dq . 



(3.28) 



W{t,x) 



The top equality is the result of Lemma 2(ii). The lower equality follows from an inte- 
gration by parts. For example 



/ g{T\lo)-Jg{T\lo)dT = 2 [ p{t\Iq) ■ q{T\lQ)dT + {q[ ■ p^^ - q\ ■ p\) 
'o Jo 
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where (^oiPo) = and {ql,p\) = k- The dr integral is that of p ■ dq over the curve 
T^{lQ,lt)- The sum of all the surface terms from ( p.27| ) is just the integral over the two 



chords, C(ro, /q) + C*(/t, r^). Thus one has that the time integral ( p. 271) is the loop integral 

The 7^2 dr term may be rewritten in loop form by employing the Poincare-Cartan 
identity in Lemma jl^. Setting the initial curve to be 70 = C(ro, ^o), the identity (|C.l 



IS 



{H{lo)-H{ro))t= I p-dq. 

J C(roM+T+ilo,lt)+git\Cilo,ro))+T- (rt,ro) 

Adding all the loop integrals gives ( ^.25| ). 

The amplitude expressions result from Mt{xQ{t,x)) = x, and the fact that the ex- 
ponentiated integral in ( |3.5| ) vanishes. This latter is a consequence of tr V^Vyi/ o / = 
tri7"(-|j) = 0, etc. □ 

The Hamilton-Jacobi equation for the p{t, x) problem is 

dtS{t, x) - H{x - yVS{t, x)) + H{x + i JV5(t, x)) = . 

All the Sj{t,x) functions, provided by ( p. 251) , are solutions of this Heisenberg H-J equa- 
tion. As t nears 0, At will have one sheet. Its phase obeys, S{0,x) = So{x). 

Lemma p!0| establishes that the amplitudes Nj{t, x) and Q:j{t, x) are finite for all x G A/t. 
When X G /C(Af) the amplitudes diverge and the representations of Theorem 1 and 2 are 
not applicable at these caustic points. If x is not in Df, it is classically forbidden. In this 
sector, U(t,x) and p{t,x) are 0{h°°). 

The area representation of S(t, x) is a consequence of Stokes' theorem. 



S{t,x) = So{xo{t,x)) + dqAdp. (3.29) 

Ja{L{t,x)) 

Above, a{L(t,x)) denotes any oriented surface with boundary L(t,x). The specific form 
of ( p.29| ) resulted form choosing the 70 curve to be C(ro,/o)- Other natural choices for 
the curve connecting vq and Iq do occur. Consider the problem where po{x) is given by 
the projection operator po = \ipo){ipo\. As was noted in Sec. 2 the initial phase function 
So{x) is the loop C{lo,ro) + 7A(ro,/o) integral (|2.13|) . By choosing 70 = 7a('"o,^o) in the 
identity (C.l) one can incorporate the initial phase into the final loop integral. So for this 



particular initial state, the entire Heisenberg phase has the form of a symplectic area 



S{t, x) = dq A dp , 

Ja{L'(t,x)) 

where L'(t, x) = g(t\-fx{ro, Iq)) + C{lt, n) . 

The phases x) and S{t,x) form representations of the X2 Lagrangian manifolds 
A(, via y = V$(t, x) and y = VS{t,x). In turn these phases define, cf. ( |2.8|) , a canonical 
transformation between the dynamical left right coordinates /(t|/o) and r(t|ro) . The At 
phases also generate dynamical flow. In the Heisenberg case one has 

/(t|/o) = 9{t\lo) = x- UVS{t, x) , r{t\ro) = g{t\ro) =x+ UVS{t, x) . 
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The presence of symmetries add structure and simplify these semiclassical expansions. 
Two Weyl symbol specific symmetries are: affine canonical covariance and fi-parity in- 
variance. 

Consider the first of these. In the notation of Appendix A, let V be the unitary 
operator defined by the affine transformation A : xi ^ Xi- Let Uyit) and py(t) be 
the symbols of the l^-similarity transformed evolutions Uv{t) = VU{t)V~'^ and pv{t) = 
Vp{t)V^^. The semiclassical covariance statements are 

These identities follow from the exact covariance property, e.g. pvif) = p(f) ° A^^, and 
the fact that asymptotic expansions have unique coefficients. 



The h parity property is specific to Heisenberg evolution. Recall that, cf. (2^), the 
symbol po{x;h) for the initial state po = \4'o){4'o\ is an even function of h. For most 
problems of physical interest this h-even property is maintained under time evolution. 
Specifically, one has 

Lemma 4. (Evolution h-parity). Suppose the real valued H{x] K) and po{x] K) have even 
h parity. Then the solution of ^l-W , p(t,x), is real and has even h parity. The associated 
semiclassical expansion has the improved convergence 

pit,x)= p''it,x)[l + Oih')], 

p'%t,x) = \det\/Mt{xoit,x))\''^/'^aoixoit,x))cos[S{t,x)/h]. 

Proof. The function p(t,x) = p{t,x;h) is defined as the unique solution of (|1.2|) . Since 
the Moyal bracket and H{x; h) are h even, the equation of motion ( |1.2| ) is invariant under 
h —>■ —h. Because initial state the po{x] h) is even it follows that p{t, x; K) is h even. □ 

It is useful to review to what extent the use of the secondary phase approach is crit- 
ical in obtaining Theorems 1 and 2. Marinov obtained the small time, $o = 0, version 
of Theorem 1 without using the embedding of Xi i^ito X2 ■ In this approach (see 
for related details) one starts with problem ( |1 . 1| ) , and assumes a WKB representation 
for U{t,x) ~ U'^'^{t,x). After applying the definition of the * product to H * U^'^{t) the 
Hamilton- Jacobi equation equation ( |3.23D for ^(t, x) is derived from the h ^ Q limit. 
The transport equations for amplitudes N{t, x) are obtained from (pTTp by extracting 
the ordinary differential equations generated by the higher order h term identities. A 
similar approach would work for the x) problem. What this strictly primary phase 
space approach misses is the recognition of the Lagrangian manifold A^ as the geomet- 
ric structure that determines the Schrodinger and Heisenberg WKB expansions. Given 
the SPS setting one can make the WKB representations global by following Maslov's 
approach. In Theorems 1 and 2 the base Lagrangian plane was always taken to be R^" 
with commutative coordinates {xa}'^- The caustic and non-diffeomorphic behavior of 
the Hi projection is avoided by replacing Hi by a different and non-singular projection: 
H*At = A*. Here A* is a a new Lagrangian plane with mixed (but commuting) coordinates 
{xa{i),yi3{i) :« = !,..., 2n}. A related point concerns the amplitude representations. Only 
in the SPS setting is it possible to express the amplitudes as determinants of a Poisson 
bracket. 
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4 Applications and Examples 



In this section we collect various interpretations, extensions and applications of Theorems 
1 and 2. For example, in both the Schrodinger and Heisenberg pictures, time evolution is 
implemented by a one parameter group. Define T(t)pQ = U{t)pQU{t)K These two group 
properties are U{ti+t2) = U{ti) U{t2) and V{ti+t2) = T{ti)oT{t2), respectively. What are 
the symplectic phase identities generated by these groups? Another question of interest 
is to establish the mutual consistency of the U'^'^{t,x) and p'"^{t,x) approximations. 

The systems ( |1 . ![ ) and ( |1.2| ) admit exact solutions if the Hamiltonians are quadratic. 
In this special case, the contrasting behavior of U (t, x) and p{t, x) is explored. Finally, 
we show how the semi-classical solutions of ( |1.2| ) behave when po{x) is not a rapidly 
oscillating function of H. 



4.1 Phase Additivity Identities 

Phase additivity identities, mirroring the group properties of U{t) and r(t), are a direct 
consequence of the loop structure of the BC problem. Denote the BC loop integrals by 

ip{t^x;L)= (b p ■ dq , A{t,x;W)= (p p ■ dq . 
JL{t,x) Jw{t,x) 

First consider the Schrodinger problem. Let ti and t2 be successive time displacements. 
For X G J\fti+t2 the BC solution determines, xo(ti + t2, x), (where the j label is suppressed). 
The loop for the full time interval is 

Li+2(ti + t2, x) = T+(/o, /1+2) + C(/i+2, ro) + C(ro, k) ■ 

Here Z1+2 = g(ti + hl^o)- As Fig. 5 shows this loop can be decomposed into the sum of 
two (or more) sub-loops defined for ti and t2, e.g. li = g{ti\lo). Let Li and L2 be the 
loops with endpoints (/q, h, ^o) and (/i, /1-1.2, tq). It follows from their definitions that these 
loops are additive, 

^1+2(^1 + t2,x) = Zi(ti,Xi) + L2{t2,x) , Xi = hli + ro) . 




Figure 5. Schrodinger loop addition. Figure 6 . Heisenberg loop addition. 

As Fig. 6 makes clear there is a corresponding loop additivity in the Heisenberg BC 
problem. Note that there is no requirement that the intermediate midpoint {ti,xi) be 
non-focal. Using the phase definitions above one obtains the following. 



17 



Proposition 3. Let x G J^ti+t2 o,f^d suppose, L\^2 = L{ + L2; = ^1 + 

are loop decompositions consistent with the BC solutions Mjj_|_t2(a;o(ti +^2,2;)) = x and 
Mt^+t2{xl{ti +t2,x)) = X, j = 1, . . . ,N, then 

ip{ti+t2,x;L{^2) = ^ih,xi;L{) + (p{t2,x;Li) , (4.1) 
A{ti + t2, X- W(^2) = Ah, x{; W() + A{t2, X- Wi) . (4.2) 

Consider the short time regime where Proposition 1 holds. Let $0 = and xo{t,x) = 
M-\x) be the BC solut ion. In this sector, it is possible to compare the rule ( [4.1|) with the 
phase addition formula established by Marinov |10[. Denote by X{t,x) the two segment 
loop T~^{xo, g{T\xo)) + C{g{T\xo),xo), where xq = lo = M~^{x). Suppose ip{T,x;X) is 
the corresponding loop phase. In terms of these two sided figures, see Fig. 7, the loop 
additivity is realized by 

^1+2(^1 + t2, X) = Xi(ti, Xi) + X2{t2, X2) + C{lo, h) + C{k, I2) + C{l2, lo) . 

The implied phase additivity now reads 

ip{ti + t2, x; X1+2) = 'f{ti,xi; Xi) + Lp{t2, X2\ X2) + ^3(2;, Xi, X2) . (4.3) 

Equation (|4.3| ) is the version of additivity found by Marinov. It is the use of the two sided 
loops X, rather than the three sided loops L, that accounts for the extra triangular area 
P3. 




Figure 7. Marinov loop addition. 

A somewhat similar additivity links the Schrodinger and Heisenberg phases with 
each other. On the left branch of Fig. 2, ip{t^xi\Xi) is defined by a two sided loop, 
Xiit^xi) = T~^{lo,lt) + C{lt,lo). Similarly, ip{t,Xr; Xr) is the right branch area defined 
by Xr{t,Xr) = T+(ro,rt) + C(ri,ro). Recall from ( |A.5| ), that the four chord loop with 
midpoints {x,Xr,xo,xi) has symplectic area P4{x,Xr,xo,xi). With this notation one has 

Proposition 4. For x E Aft let (x, x^, x^, xj) be the midpoints of the Heisenberg BC loop 
solution {W^(t,x)} I . Then 

A{t, x; W^) = ipit, xi; xi) - ip{t, xi; Xi) + P^{x, x^ xl x\) , j = 1, ■ ■ ■ , iV . (4.4) 

Proof. This follows directly from the loop decomposition 

W{t,x) = C(ro,/o) + r+(/o,/t) + C(/i,n)+T-(n,ro), 

= X(ti, Xl) - X{t2, Xr) + C(ro, /o) + C(/o, k) + C{lu rf) + C(n, h) . □ 
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4.2 Mutual Consistency of U'%t) and p'%t) 

One may build the Heisenberg evolution p{t, x) from its Schrodinger components, via 
p{t,x) = U{t) * po * U{t). The goal of this subsection is to show that the use of the 
approximation U{t, x) ^ U^^{t, x) combined with the evaluation of the * product by the 
stationary phase approximation constructs p'^'^{t,x). This calculation will be carried out 
using the single sheet BC solutions found in Proposition 1. All the amplitudes and phases 
of U^^it^x) and p'^'^it^x) will be written without the multi-sheet subscript, j. This type 
of composition argument was used by Berry and Balazs to construct p^'^(t,x) in the 
case where is one dimensional. 

First observe that the required form of the is the Marinov version obeying 

the initial condition U{t,x) = 1. Since $o = 0, the amplitude takes the simplified form 

N{t,x) = [det'^{l + Vg{t\Mi\x)))]-'^\ (4.5) 
Employing the three function * product ( |A.8D , the approximate p{t,x) is 



[U'''{t)*po*U^'{t)]{x)=l3{t,x) 

= cj U'\t, xi)p{x2)U'^{t, xa) 5{xi -X2 + x^-x) e^^^^^^'^^'^^^/^dxidxzrfxa , 

= c ! N{t, xi) N{t, X3) aoix2) e*M*'^i)-^(*'^3)+So(^2)+P3(xi,x2,x3)]M ^^^^^^ _ (4^5) 



where, X2 = Xi + X3 — x and c = {nh) 
Denote the total phase in ( [4.6|) by 



0(xi, X3; t, x) = ip(t, xi) - (fit, X3) + So{xi +X3-X) + Psixs, X, xi) . (4.7) 

Let ^ = {xi,X3) and Q"{^]t,x) be the Hessian of 9 with respect to the ^ variables. If 
xi = xi{t,x) and 0:3 = 0:3(^,0;) is an isolated non-degenerate critical point of 6, i.e. the 
solution of V59(^; t, x) = 0, then the standard stationary phase approximation to (|4.6|) 
states 



if'it^x) = 22"Ar(t,xi)Ar(t,X3)ao(x2)|det0"(xi,a;3;t,x)r^/2 



X exp 



1 

-Q{xi, X3; t, x) + — sgne"(a;i, X3; t, x) 
n 4 



(4i 



Above, the functional form of the constraint condition is X2(t, x) = Xi{t, x) + X3(t, x) — x. 
In component form, the critical point equation reads 

xi = X + iJVv2(t,a;3) - |JVS'o(xi +X3 - x) , (4.9) 
X3 = x + iJV(/?(t,Xi) + |JVS'o(xi+X3-x). (4.10) 

For small times and under the same assumptions as Proposition 1, one can show that 
the system above defines a contraction mapping and so has unique solution xi = Xi{t, x), 
= X3(t,x). In fact the set of points {x2(t,x),xi(t,x),x,X3(t,x)} are exactly the four 
midpoints {xo,xi,x,Xr} appearing in Fig. 2. The constraint condition for X2 permits one 
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to replace the P3{x3,x,xi) area by the P4{x3,x,xi,X2), cf. ( [A .91 ). Substituting these 
critical point arguments into the total phase and using Proposition 4 one has 



Q{xi(t, x),X3(t, x); t, x) = S(t, x) . 

This demonstrates that the stationary phase produced by U^'^{t) * pq* U^'^{t) is the same 
as the Heisenberg WKB phase S{t,x). 

Verifying amplitude consistency is more elaborate. A non-trivial determinant identity 
for O" in terms of the individual left-right flows is needed. Let Ki = Vg{t\l), I = X2{t, x) — 
|V JS'o(x2(t, x)) and Kr = Vg(t\r),r = X2(t,x) + |VJS'o(x2(t, x)). Using identities 
obtained from (|4.9| ) and ([4.10|) one can establish 

det 9"(xi(t, x), Xsit, x); t, x) 

= 26"( det(l + Ki){l + Kr))~' det [Ki{l - lJS'^{x2)) + Ml + lS'^{x2))] . (4.11) 

Putting ( [4.11|) into (|4.^ ) recovers the Theorem 2 amplitude a{t,x). In the small time 
sector all the determinants above are positive and sgnQ"{xi, X3;t,x) = 0. Altogether one 
has the identity 

[W'{t)*po*U^fP\t,x)=p^%t,x). (4.12) 



4.3 Quadratic Exactness 

If the Hamiltonian is a quadratic function, then the semiclassical quantities U^'^{t,x) and 
p^'^{t,x) are exact solutions of ( |1 . 1| ) and ( |1.2| ). Thus it is of interest to obtain the explicit 
formulas for the phases and amplitudes which arise is this special case. 

The general time dependent Weyl symbol Hamiltonian of quadratic form is 

H{t, x) = ^x- H"{t) X + H'{t) ■ X , 

where H"{t) is a 2n x 2n symmetric matrix and H'{t) is a 2n component vector. Both 
H"{t) and H'{t) are real and t-continuous. The quantum Hamiltonian, H{t), results if x 
is replaced by x in H(t,x). 

The classical trajectories g{t, s\x) generated by a quadratic H{t, x) are linear functions 
of X. Explicit representations of g{t, s\x) are obtained as follows. 

The dynamical equations for g{t,s\x) and Vg{t,s\x) are the linear ODE's 

^g(t, six) = JH"{t) git, s\x) + JH'it) , (4.13) 
at 

^Vg{t, s\x) = JH"{t) Vgit, s\x) , (4.14) 
at 

with initial conditions g{s, s\x) = x and Vg{s, s\x) = I. Both ( [4.13| ) and ( |4.14| ) are Jacobi 

field equations, but ( [4.14| ) is homogeneous, while ( [4.13|) has an inhomogeneity, JH'{t) . 
Let K(t, s) be the Sp{2n) solution of 

^K{t,s) = JH"{t)K{t,s) , K{s,s) = I. (4.15) 
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The solution of the system ( |4.15| ) is unique, which imphes that Vg{t,s\x) 
terms of K{t, s), the forward {t > s) and backward {s < t) flows are 

g{t,s\x) =K{t,s)[x + F{t,s)], 
g{s,t\x) = K{s,t) X — F(t, s) , 



F{t,s) = I dTJK{T,sfH'{T) t>s. 

J s 



K{t,s). In 

(4.16) 
(4.17) 

(4.18) 



The solution ( ^.16] ) is verified by putting it into ( |4.13| ) and using the symplectic identity 
K{t,s)JK{t,sf = J. From ( ^l6D one obtains j ^Tft) by using the composition laws 
g{t, s) o g[s, t) = Id and K{s, t)K(t, s) = L 

The behavior of the Heisenberg evolution is summarized in the following lemma. Note 
that for time dependent Hamiltonians the boundary condition problem needs to be stated 
in terms of g{t, s). So the left side of ( |3.14| ) has two times, e.g. Mt^s{x'). 



Lemma 5. For quadratic Hamiltonians the semiclassical density matrix is the exact so- 
lution of p'^^{t,x) = p(t,x). The structure of p'^'^(t,x) is given by: 

(i) The BC problem Mt^s{xo{t,x)) = x has a unique global solution given by the back- 
ward flow 



Xo{t, x) = K{s, t) X — F{t, s) 



(4.19) 



S{t,x) = Soixo{t,x)) 



det VM 



t,s 



X 



(iii) 



p"'{t,x) = po{xo{t,x)) e 



iSo(xo(t,x))/h 



(4.20) 



(4.21) 



Proof . The x-linearity of ( |4.16|) combined with ( 3.14 ) implies ,,(x') = g{t,s\x'). Since 
g{t,s) has inverse g{s,t), ( [4.19|) holds. The phase identity in ( [4.20|) follows from the fact 
that chords evolve into chords, specifically g{t, s\C{lo,rQ)) = C{lt,rt). For this reason the 
Heisenberg loop L(t,x) collapses to zero, so that all the WKB phase resides in the 5*0 
term. The determinant identity is a consequence of VMt^s {x') being a symplectic matrix. 
That p'^^it.x) = p(t,x) is verified by direct substitution into (|1.2|) . □ 
The X, y representation of the manifold At is 

y{t, x) = VSo{K{s, t)x- F{t, s))K{s, t) . 

This formula shows that is single sheeted. Statement (ii) establishes that is non- 
singular for all t. 

Next consider Schrodinger evolution. Let a{s) = {t G ]R| det(i^'(t, s) — I) = 0} denote 
the set of caustic times where the amplitude function of U^^{t,x) diverges. The behavior 
of Schrodinger evolution is summarized below. These are known results ([|l( 
their proof parallels that of Lemma |^. 
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and 
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Lemma 6. For quadratic Hamiltonian Schrddinger systems with initial data U{0, x) = 1, 
the semiclassical solution of is exact at all non-caustic times. The structure of 

U'^^{t,x) = U(t,x) fort G ]R/cr(s) is given by 

(i) The BC problem Mt^s{xQ{t,x)) = x has a solution given by: 

xo{t, x) = {Kit, s) + I)-\2x - Fit, s)) , (4.22) 

(ii) 

iV(t,xo) = 2"|det[ir(t,s) + J]|-^/^ (4.23) 

(iii) 

$(t, x)^-x. . - Fit, s) X + at, s) , (4.24) 

where C(t, s) is an x-independent constant. 

The manifold is determined by the identity 

yit, x) = -2 J-^^^iy^ [iKit, s)-I)x + Kit, s)F(t, s)] . 

The x-hnearity of yit, x) means that Aj is a hyperplane. For each t^. G cr(s), where Kit^, s) 
has a -1 eigenvalue, A^^ is singular and | yitc,x)\ = oo. As t ^ t^, U'^'^it,x) has a delta 
function behavior. 

In comparing these Heisenberg and Schrodinger solutions, it is evident that the main 
difference is that pit, x) is always caustic free, while Uit, x) has singular behavior at t = tc- 



4.4 Evolution of semiclassically admissible observables 

An important special case of the Heisenberg evolution occurs when the initial operator 
Po is semi-classically admissible, and as a result Sq = 0. Most observables of physical 
interest (energy, momentum, angular momentum and their products) are semi-classically 
admissible. 

How does Theorem 2 simplify in this case? First note that Sq = implies that the 
initial Lagrangian manifold Aq = ffi^". Furthermore, the BC equation ( pj.l| ) becomes 
Mf(x') = git\x) = X, and so Mf^ix) = gi—t\x) for all t,x. From this it follows that 
det VMt(x') = 1 and so there are no caustics. Since Sq = 0, the left, right and midpoint 
BC solutions are the same Xoit,x) = loit,x) = ro(t,x). This means the left, right, and 
midpoint flows collapse to a common trajectory and as a result the Heisenberg loop L(t, x) 
reduces to a point. As a consequence S'(t,x) = 0. Altogether one has 

p'%t,x)=poigi~t\x)). (4.25) 



This simple transport expression for p^'^it, x) is the point of departure |T^, |T3|, [Tj, |T5|, IT 



for higher order h expansions of p(t,x). Identity (|]2|) is an example of Egorov's theorem 



The fact that Sit,x) remains zero is in marked contrast to the behavior of U^^it,x) 
where the phase $(t, x) spontaneously becomes non-zero even if its t = value van- 
ishes. This Heisenberg null phase stability is a consequence of the symmetry 7^2(2;, y) = 
—7^2(3^, —y)- In turn, this symmetry means that yit) = is a, constant of motion for X2 
trajectories that have initial value (xo,yo) = (^^o^O). Equivalently, A^ = Aq. 
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5 Conclusions 



The representations of operators by Weyl symbols gives one a complete statement of quan- 
tum mechanics set in phase space. The Wigner transform is a Lie algebra isomorphism 
between operators (with commutator bracket) and symbols (with Moyal bracket). In this 
formalism, the Schrodinger and Heisenberg evolutions, U{t,x) and p{t,x), are generated 
by Hamiltonians that are pseudodifferential operators. These two problems are made 
amenable to semiclassical approximation by embedding them in an extended phase space 
X2 = T*xi- Set in X2, the two evolution problems become special cases of an extended 
Schrodinger evolution. 

The geometry of the WKB approximation is controlled by At, the dynamical La- 
grangian manifold in X2- The associated BC problem is solved in terms of loops built 
from chords and trajectories in PPS. The symplectic area of these loops form the WKB 
phase functions. The amplitudes are X2 Poisson brackets constructed from the family of 
constraints that define the initial and final BC Lagrangian manifolds, Aq and N't C M^"-. 

An alternate characterization of Heisenberg evolution in symbol space is found in 
Marinov's path integral representation of p{t, x). This path integral has the following 
structure. Write the mapping po{x) ^ p{t,x) in terms of an integral transform 



p(t,x) = J dx' K{t;x,x') po{x') , (5.1) 
where the kernel K is /i~"Tr[t/(t) A(a;')?7(t)^A(x)]. Define an action functional by 
S{t;x,y)= [ dT[y{T)-S:{T)-n2{y{T),x{T))]. 



The path functions x{t) obey the boundary conditions, x(0) = x' and x{t) = x, while the 
functions y(r) are unrestricted. Then K is the path integral. 



(5.2) 



Here Dx Dy are infinite products of the dimensionless measures h~"'dx and h'^^dy. 

A second perspective of ( |5.2|) follows from the extended Schrodinger equation ( p.l|) . 
In the SPS setting, p{t, x) may be interpreted as a wave function generated by the 7^2 
evolution in L^(i?^"). The kernel K is then the xi coordinate space Dirac matrix element 
(a:|exp (^—i'H2t/h)\x') and representation ( |5.2| ) is the corresponding Feynman path inte- 



gral. The parallels with the p^'^{t,x) approximation are evident. The functional S{t;x,y) 
has the same form as the generic WKB phase (|3.3| ), but with ^(r), x{t) functional argu- 
ments replacing the classical flow G{t\ x',y'). 

The symplectic area WKB structure found for the evolutions U{t) and p{t) seems 
to be universal. An alternate version of Theorem 1, suitable for quantum systems with 



external electromagnetic fields, has been recently published in . That paper introduced 
a gauge invariant * product of the Berezin type. The {xi,X2,X3) triangle phase of this 
magnetic * product is a symplectic area with respect to the Faraday 2-form. The short 



time electromagnetic generalization of our Theorem 1 is found in Theorem 6 of p2] . Using 
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the SPS methods developed here combined with the with symbol calculus of the magnetic 
* product, the extension of Heisenberg WKB representation of Theorem 2 to a gauge 
invariant form is straightforward. 

Another direction of generalization is to consider non-Weyl g,p quantum orderings. 
The simplest of these alternate ordering schemes [§] is 

1 3 2 

fiQ,p) = /(a Q +(1 - a) q,P) ■ 

where a G [0, 1]. Normal, anti-normal, and Weyl ordering occur for parameter values a = 
0, 1, |. For all orderings with a G (0, 1) there will be symplectic area WKB representation 
of Schrodinger and Heisenberg evolutions. 

The extended phase space method developed here is expected to be a good platform 
for further study. Since the X2 Schrodinger evolution ( |2.1D is unitary in L^(M^"'), norm 
error bounds for \\p(t) — p^'^{t)\\2 may be obtained with standard 0, ^ techniques. Higher 
order h expansions, with U^'^{t, x) and p*^(t, x) as first terms, are straightforward to derive 
using the known higher order transport equations. Similarly the available globalization 
procedures 0, ^ readily apply in the X2 framework. 



A The Star Product 

Comprehensive overviews of the Weyl symbol calculus are found in the books [18, 33]. In 
this appendix we summarize just the features of pseudodifferential operators and the star 
product that our semiclassical analysis requires. It is useful to begin with Weyl's original 
idea of quantization Let / be the fi- Fourier transform of / : T*M'^ —>■ C 

fix) = [ dy f{y) exp (^y ■ x) ^ {Fn,y^J) [x) , (A.l) 

then the Weyl quantization of / is 

/ = / dyf{y)exp(^y-x)=f{q,p). (A.2) 

In this approach the family of unitary operators exp ^y ■ x is associated with the exponen- 
tial functions exp j^y ■ x. By Fourier superposition this association is extended to construct 

a general operator /. The notation / = f{q,p) specifies a Weyl ordered function of the 
canonical operators {q,p). 

Upon first inspection the quantization procedure above does not clarify if there is a 
well defined inverse: given an arbitrary operator / on 7i how is f{x) determined? Moyal 
solved this problem by showing that the inverse map / i-^ / is given by the Wigner 
transform 

fix) = IfUx) ^ [ dve-^-^'\q + \v\f\q - \v) . (A.3) 



Although there are a number of other consistent schemes |14] which associate functions 



on phase space with operators, normal ordering [p acts before g]. Wick ordering [{q + ip\ 
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acts before {q — ip)], etc., Weyl quantization has the advantage that 1) it treats q and 
p symmetrically, 2) self-adjoint operators have real symbols, and 3) the Moyal bracket 
is an even function of h, in particular the semiclassical correction to the leading Poisson 
bracket term is 0(^^), not 0{h). 

A basic symmetry of the Wigner-Weyl isomorphism is its covariance with respect to 
affine canonical transformations. These transformations are parameterized by a matrix 
R G Sp{2n) and a displacement xq , 

X = A(x) = R~^{x — Xq) , A~^{x') = Rx' + x'q . 

The pair (-R, Xq) also defines a family of unitary operators, V = V{R,Xo) obeying 

VxV-^ = A{x) = R-^ {x -xqI). 

If T{xq) is the Heisenberg translation operator {T{xq)~^xT{xo) = x + XqI) and M{R) 
is a metaplectic operator (i.e. a unitary operator obeying M~^xM = Rx) then V = 
T{xq)M{R). Affine covariance relates the Weyl symbol of a similarity transformed 
operator to its A transformed symbol. Using the density matrix as an example, this 
universal covariance property |T5[ is 



{Vmv-')^ = = P(t)w o A~' . (A.4) 

In classical mechanics, an observable f{x) is a real valued function and the product 
of observables is the commutative product of functions /i(x)/2(x). In the Weyl symbol 
version of quantum mechanics, the commutative product becomes the non-commutative 
star product. There are three distinct types of representation of this star product; 1) 
the integral definition, 2) the left-right operator characterization, and 3) Groenewold's 
exponential derivative formula |TB|, |3^. Each of these three forms is useful in a different 



way. The integral form is convenient for proving the associativity of the * product and 
determining the form of multiple products; the left-right forms represent the * composition 
as a pseudo-differential operator (Lemma ||), and the Groenewold formula gives explicit 
coefficients for the small /i-asymptotic expansion of the star product. 

Consider first the integral representation. Let fj, j = 1, . . . , N he N different opera- 
tors. The Weyl symbol of their product is denoted by 

[/1/2 ■ ■ ■ /7v]w = fi* f2* ■ ■ ■ * In ■ 
The higher order * product requires the following two functions 

Af-l N 

Pn{xi,X2, . . . , xn) = 2 ''^^(—lY'^^'^^Xi A Xj , Xi A xj = Xi ■ Jxj , (A. 5) 

i=l j>i 

SNixi,X2, . . . , Xn) =Xi-X2-\ h (-l)^+^XAr . (A. 6) 

Lemma 7. Let fj G L-^{M.'^^) be Weyl symbols, j = 1,...,N. For even N > 2, the 
integral form of the * product is 

(/i * /2 * ■ ■ ■ * fN)ixN+i) = c^/^ j dxi... dx^fiixi) ■ ■ ■ fwixN) exp [iPN+i{xi, XN+i)/h] ■ 

(A.7) 
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where c = (nh) . For odd N > 3, 

(/l * /2 * ■ ■ ■ * fN){xN+l) = C^^"^)/^ j dXi... dXNf l{Xi) ■ ■ ■ fn{xN) 

6 (S'Ar+i(xi,...,X7v+i))exp [iPN{xi,. . . ,XN)/h\. (A. 8) 



Proof. The result for N = 2 follows directly from (|A.1| ) and (|A.2|) . The higher order 
products follow from the recurrence relation 



Pn{xi, . . .,xn) = Pn~i{xi, . . .,xn-i) + (-1) ^ 2xn A Sn-i{xi, . . .,xn-i) ■ (A. 9) 

and an induction argument. □ 

The phase Pat in the * product has a simple geometric meaning. It is the action 
line-integral 



P/v(xi, . . . , Xat) = (p p ■ dq = dqAdp. (A.IO) 

The geometric figure in the integral is an A^-sided polygon in T*W^ having (xi, . . . ,xn) as 
the midpoints of its successive sides. Note that in Lemma the phase P/v is needed 
only for odd. When A^ is odd, this polygon is uniquely determined by the mid- 
points (xi, . . . ,xn)- Denote the boundary of this polygon by with orientation xn 
xn-1 ■ ■ ■ — > xi — Xn', Sat is any oriented surface with outer boundary L^. Thus P/v is 
the symplectic area of this polygon. For a single star product fi * f2 the polygon is a 
triangle and P3 is its symplectic area. 

A variety of symmetries for the symplectic area, P/v, follow from its definition ( |A.5| ). 



PNixi + a, . . . ,xn + a) = Piv(xi, . . . ,XAr) - 5iv2a A S'Ar(xi, . . . , Xiv) , (A.ll) 
Pn{-xi,-X2, . . . ,-xn) = Pn{xi,X2, . . . ,Xn) , (A.12) 

Pn{x2,Xs,. . . ,Xn,Xi) = PAr(xi,X2, . . . ,XAr) , (A.13) 

Pn{xn,xn-i, ■ ■ ■ ,xi) = -Pn{xi,X2, . . . ,xn) ■ (A.14) 

The above identities express translation, refiection, cyclic and anti-cyclic permutation 
invariance. The quantity 6n is 1 if A^ is even and zero otherwise. 



These geometrical results for A^ = 3 are due to Berezin |^ . It is now known [g^ that 
the * product on all symmetric symplectic spaces has a phase that is a triangle related 
symplectic area. Other representations of the A^-order * product are found in references 
0,11 

A second perspective on the * product is to understand it as a \I'DO. These operators 

are generated by smooth SPS functions, which are also called symbols. A \1/D0 acting on 

2 1 

the function f(x) in denoted by T-C{X,Y)f{x). This operator is normal ordered, which 

means that Y = —ihVx acts first followed by X, which is multiplication by x. The 
2 1 

operator 'H{X, Y) is formally defined as follows, 

[n{X, Y)f] i^) = [FH,y^,n{x, y)Fn,,^yf] (x) . (A.15) 
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2 1 

The presence of h in the Fourier transform Fn^^^y means that Ti.{X, Y) is a function of h, 
and that the h ^ asymptotics is easy to study. Passing from a normal ordered symbol 
to a Weyl ordered symbol is easy and is achieved by the transformation 

T-(.Weyi{x, y) = exp(^V^. ■ Vy)n{x, y) . (A.16) 

In this paper the symbols with SPS arguments (x, y) generally correspond to normal 
ordered operators, while symbols with the PPS argument x correspond to Weyl ordering. 

2 1 

The definition of H{X, Y) is completed by placing H in a suitable class of functions. 
Two important function classes enter in Maslov's proof of ( ^.11 ). The first is T'^{M.f^). 
This class denotes those C°°{M.f^) functions Ti which satisfy the uniform growth estimates 

\VlVP^n{x,y)\<C,,{l + \x\r{l + \y\r, |7|,|p|>0, (x, G x Mf = , 



where integer m > and the positive constants C^p are independent of {x,y). The 
symbols in the class T™(]R^") are h independent. Now widen this class so as to in- 
clude semiclassically admissible operators. These operators depend on h and admit an 
asymptotic expansion about h = 0. One says H{x,y; h) belongs to the class T^(M^"), if 
n E C°°(Mf X R+) and has the h expansion 

J 

n{x,y; h) = Y,mm,{x,y) + h'+'nj{x,y; h) . (A.17) 

j=0 

The coefficient functions Tij G T™, and the remainder term TZj has the {x, y; K) uniform 
bound 



d 



\VIVI{—Ynj{x, y- h)\ < + |a;|r(l + \y\r , ItUpI, A > . 

Now return to the * product. Its left-right multiplication forms are the following. 
Lemma 8. For all /i, /2 G 5°°(R2") 

/i*/2(x) = /i(l-|JF)/2(x), (A.18) 
/i*/2(x) = f2{X+\JY)h{^)- (A.19) 
Proof. Make a change of variables in ( [A. 7] ) so that 



/i * f2{.x) = {2-nK) 2" j dx2 dy2 fi{x - \ Jy2) f2{x2) exp (^^y2 ■ {x - X2, 

2 1 

In the integral one can replace fi{x — l/2Jy2) with fi{X — 1/2J Y)- Moving this last 
factor outside the integral gives ( [A.18| ). A similar argument gives ( |A.19| ). □ 



Further review of the left right operator method in quantization and the use of the 
left right projections to provide representations of SPS is found in [ pT^ . 
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B Jacobi Field Applications 



The Jacobi field of a classical system characterizes its stability. Given a trajectory g{t\x), 
the 2n x 2n matrix solution of the linear system 

f^^gm - JH'\g{t\x))Vg{t\x) = , (B.l) 
Vg{0\x) = I, (B.2) 



defines the associated Jacobi field. The solutions of (|B.1| , p.2|) are symplectic matrices. 
For small times Vg{t\x) will not differ greatly from its initial value I. This is made explicit 
in the following estimates. The norm ||-|| denotes the operator norm on M^". 

Lemma 9. Let the Hessian of H have the x-uniform hound \\H"{x)\\ < Ci < oo, then the 
Jacobi field has the growth estimates 

\\Vg{t\x)\\ < e^^* , \\Vg{t\x) - I\\ < (e^^* - 1) . (B.3) 

Proof. Define ^{t\x) = ||V^(t|a;)|| and use |||V^(t|2;)|| < |||V^(t|a;)||. Thus equation 
( p.l|) and the Hessian bound implies 

< \\JH"{g{t\x))Uit\x)<c,at\x), 
Now estimate ||V(7(t|a;) — /||. The integral equation equivalent of ( p.l| , |B.2]) is 



Vg{t\x)=I+ [ JH"{g{T\x))Vg{T\x)dT. 
Jo 

Taking the norm of this relation gives 

||V^(t|x) - J|| < [ \\JH"{g{T\x))U{T\x)dT <ci [ e"'* dr = (e^^* - 1) . □ 
Jo Jo 

Next we consider the detailed linkage between the BC problem, finite WKB ampli- 
tudes, and the non-singular set of with respect to the projection Hi. Although this 
connection is implicit in the Maslov-Fedoriuk book 0, its precise statement is important 
for the BC solvability analysis and the interpretation of the WKB amplitude singularities. 

First let us collect several definitions and facts. Recall that Dq is the support of Sq 
or $0 • Furthermore the projection Hi : Aq — * Dq is diffeomorphic, with explicit inverse 
n^i(x') = {x',VSq{x')) , x' e Dq. The regular subset of A* = G{t)Ko is Tit = At/S(At). 
Let Tm{At) denote the tangent plane at m G A^. The statement, m e TZt, means that 
no vector in Tm(A() has a zero Hi projection, i.e. dim IliTm{At) = 2n. Since A^ is a 
Lagrangian manifold, Tm{At) is a 2n dimensional symplectic vector space for each m. 
With this notation we have 
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Lemma 10. Let x' E Do = HiAq , then 

det VMt{x') ^0 ^ G{t\x\ VSq{x')) e Tit 



(B.4) 



Proof. The X2 Jacobi field VG{t\m') maps vectors on Tm/(Ao) into vectors on TG{t\m'){,At) ■ 
The matrix VG(t|m') is symplectic, specifically det VG(r|m') = 1. Let {cj}^" be a basis 
of Do ^ M^n ^iti^ origin at x' . Then {W^i = (ci, S'^'(a;')ei)}?'^ is a basis of T^'(Ao), m' = 
(x', V^olx')). Since VG{t\m') is invertible {VG(t|m')Vri}f is a basis of TG(t|„')(At). The 
Jacobi field for 7-^2 fiow is 



VG{T\x\y' 



-,[Vg{t\V)+Vg{t\r') 



-Vg{t\V) + Vg{t\r')]J 



J[Vgm - Vg{t\T') ] -\J[Wgm + Vg{t\T') ]J 



Consider the conditions which ensure that m G TZt- Let Z = {Zx,Zy) represent an 
arbitrary non-zero vector in Tm(At). Since {'VG{T\m')Wi}i"' is a basis one has 



2n 2n / 

Z = Y^ \iVG{T\m')Wi = J2 A*VG(r|m') ( 
j=i j=i ^ 



S'q (x )ej 



where at least one of the constants Aj are non-zero. Project Z onto M.^^., giving 
U,Z = Z, = l^ [g{t\x' - i JV5o(x')) + git\x' + i JV5o(a;'))] K , 



2n 



i=l 



Thus a vector in Tm(At) can have zero Hi projection iff det VMi(x') = 0. □ 

With the replacement of Sq by $o and Mt by Mj, the equivalence ( p.4| ) holds for the 
Schrodinger problem. 



C A Poincare— Cartan Identity 

The following identity is required in the proof of Theorem 2. 

Lemma 11. Let 70 be an arbitrary smooth curve from Xi to x^. Let (T(t|7o) be the surface 
defined by the images of under the flow (?(r, to) , tQ <t <t. Then 

/ dqAdp= (h p.dq= [H{t, g{T,to\xi)) ~ H{t, g{T,to\x2))]dT . (C.l) 

Ja{t\'yo) J da(t\"io) J to 

Proof sketch. Decompose the curve 70 into segments that are singularity free with 
respect to some Lagrangian plane. For notational convenience denote this plane by W^. 
Now locally embed this curve into a non-singular Lagrangian manifold Ao having phase 
So(q'). Denote by Xt = g{t, to|Ao) the time evolution of Ao. Extend the original phase space 
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Xi by adding time to the g-coordinate manifold, {q,t]p,pt) E Xi = T*R'^^^. Let s{t,q) 
be the phase of A^, then the manifold 

A+ = {{q,t-p,pt) e xt\{P,Pt) = {Vs{t,q),dts{t,q))} 

is Lagrangian. On A+ the 1-form uj^ = pdq — H{t,q,p) dt is closed by virtue of the 
H-J equation. Closed loop integrals on A^ vanish and can be arranged to give ( |(J.1| ). □ 
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